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ABSTRACT 



The effects of Hall electromotive forces (HEMFs) on the linear stability of 
protostellar disks are examined. Earlier work on this topic focused on axial 
field and perturbation wavenumbers. Here we treat the problem more generally. 
Both axisymmetric and nonaxisymmetric cases are investigated. Though seldom 
explicitly included in calculations, HEMFs appear to be important whenever 
Ohmic dissipation is. They allow for the appearance of electron whistler waves, 
and since these have right-handed polarization, a helicity factor is introduced 
into the stability problem. This factor is the product of the components of 
the angular velocity and magnetic field along the perturbation wavenumber, 
and it is destabilizing when negative. Unless the field and angular velocity are 
exactly aligned, it is always possible to find destabilizing wavenumbers. HEMFs 
can destabilize any differential rotation law, even those with angular velocity 
increasing outward. Regardless of the sign of the angular velocity gradient, the 
maximum growth rate is always given in magnitude by the local Oort A value 
of the disk, as in the standard magnetorotational instability. The role of Hall 
EMFs may prove crucial to understanding how turbulence is maintained in the 
"low state" of eruptive disk systems. 

Subject headings: accretion, accretion disks — 
magnetohydrodynamics — instabilities — turbulence 
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1. Introduction 

The stability of differentially rotating gas disks depends sensitively upon whether 
or not a magnetic field is present. A weak (subthermal) magnetic field undermines the 
stabilizing influence of Coriolis forces (Balbus & Hawley 1991), and the resulting turbulence 
greatly enhances internal angular momentum transport (Hawley, Gammie, & Balbus 1995; 
Balbus & Hawley 1998). This important behavior, which we refer to here as the standard 
magnetorotational instability (MRI), probably is the underlying cause of "anomalous 
viscosity" in accretion disks. While other modes of transport (e.g. waves) are certainly 
possible, the classical enhanced turbulent transport associated with a disks is almost surely 
MHD turbulence (Balbus & Papaloizou 1999). 

Protostellar disks remain problematic. They seem too dense and too cool to be 
magnetically well-coupled over their full radial extent. Understanding how protostellar and 
protoplanetary disks behave is an ongoing theoretical challenge. It has also been recently 
noted (Gammie & Menou 1998, Menou 2000) that the outer regions of dwarf novae disks 
may have a very low ionization fractions. That protostellar disks and dwarf novae both 
show episodic eruptions may be no coincidence; in both systems the question of magnetic 
coupling (and therefore the presence of MHD turbulence) is a delicate one. This paper 
addresses a very small but fundamental part of this problem: the linear effects of Hall 
electromotive forces (HEMFs) on disk dynamics. While the results are more broadly 
applicable, protostellar disks will be the focus of this paper. 

The importance of the Hall effect in modifying the MRI was pointed out in a recent 
paper by Wardle (1999; hereafter W99), who analyzed the linear stability of a protostellar 
disk in the presence of a vertical magnetic field. Such disks are generally in the Hall regime. 
That is, their combination of densities, temperatures and ionization fractions places them 
in a realm of parameter space in which departures from ideal MHD are important, which 
are strongly influenced by HEMFs. The Hall effect tends to be important simultaneously 
with finite conductivity. Numerical simulations of low ionization disks (e.g., Fleming, Stone, 
& Hawley 2000) have not yet included Hall terms, but have focused rather exclusively on 
finite conductivity. 

W99 found that standard, locally unstable MRI modes may either be stabilized or 
destabilized depending upon whether f2 • B is positive or negative. (Here, as usual, S7 is 
the local angular velocity vector of the disk, and B is the magnetic field.) In other words, 
HEMFs appear to introduce a helicity effect into the problem. 

This immediately raises an interesting question: what happens if f2 • B — 0? Would 
the Hall effect disappear? Since the W99 analysis was restricted to axial fields and axial 
wavenumbers, it could not address this problem. Clearly HEMFs do not simply vanish when 
the magnetic field lies in the disk plane, and the identification of Q • B as the fundamental 
Hall stability parameter therefore seems suspicious. Identifying this key parameter more 
precisely is important, for it is at the heart of MHD stability in low ionization disks. 

We are motivated, therefore, to examine the stability of Hall regime accretion under 
very general circumstances. We restrict neither the local field geometry, nor the local 
wavenumber direction. Furthermore, we have chosen to couch the analysis in a dynamical 
context, instead of ascribing the inductive Hall terms to effective conductivities (W99). The 
effective conductivity approach offers terseness and mathematical economy; the dynamical 
approach makes the physics of the coupling between the magnetic field and the disk 
coupling more transparent. Both methods must of course ultimately be equivalent. It is 
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enlightening, however, to view this important problem from a different physical perspective, 
and we have accordingly chosen this route. 

The plan of the paper is as follows. In §2, we discuss the basic formulation of the 
problem. This includes the various physical and compositional parameters, the fundamental 
equations, and an order of magnitude estimate for the relative sizes of the Ohmic, Hall, 
and ambipolar diffusion terms in the induction equation. Section 3 presents the linear 
dispersion relation for a uniformly rotating disk, because this is the simplest setting to see 
the interplay between rotation and the Hall EMFs. Section 4 is a detailed and extensive 
analysis of differentially rotating disks for arbitrary field geometries and perturbation 
wavenumbers. Section 5 summarizes the findings of this paper. 



2. Basics 



2.1. Physical and Compositional Parameters 

A typical protostellar disk consists of molecular gas with cosmic abundances, and spans 
a temperature range from 10 K in its outer regions to a few 10 3 K near the central star. 
The characteristic disk size may be up to of order 100 AU, and its total mass will typically 
be less than roughly 0.1 of the central star's mass. If 0.01 M & is spread over a cylinder of 
radius 1 AU and height 0.1 AU, this would yield an average mass density p ~ 2 x 10~ 10 g 
cm~ 3 , and an integrated column density of £ ~ 3 x 10 2 g cm~ 2 , fiducial numbers one should 
bear in mind. 

The dynamics of the disk is likely to be dominated by the Keplerian orbital velocity 

vk- 

2 GM 2 2 

v K = -ft- = Ru (i) 

where M is the central mass, R the cylindrical radius, and Q the angular velocity. The 
disk is assumed to be thin; the characteristic temperature T is much less than any virial 
temperature. The scale height h (defined below) satisfies h/R <C 1. 

The number density of species X is denoted n x . The H 2 and He number densities are 
assumed to be related by 

n H e = 0.2n H2 , (2) 
which implies a total neutral number density 

n = n H2 + n He = l.2n H2 , (3) 

and a neutral mass density 

p = 2.8n H , 2 m p , (4) 
where m p is the proton mass. This gives a mean mass per particle of 

p = p/n = 2.33m p , (5) 

and isothermal and adiabatic sound speeds respectively of 

4 = 0A29kT/m p , a 2 = 0.6kT/m p , (6) 
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where k is the Boltzmann constant. For simplicity we will define the disk scale height h by 
h = C5/Q. Under conditions of interest, the dominant ion will generally be once-ionized 
potassium K + , and accordingly we take a mean ion mass of 

mi = 39m p (7) 

In molecular cloud studies, rrii = 30m p is more typically used (e.g., Draine, Roberge, & 
Dalgarno 1983). 

The coupling between ions and neutrals depends upon the rate coefficient, denoted 
(av) iN . In essence, this is the the product of a cross section with the component of the 
relative velocity along the drift axis, averaged over the electron Maxwellian distribution 
(Spizter 1978). Similarly, there is an electron-neutral coupling rate, {<rv) e ]y. We follow 
Draine et al. (1983) and Blaes & Balbus (1994), adopting the values 

(av) iN = 1.9 x l(T 9 cm 3 s" 1 (av) eN = 1(T 15 ( ^K) ^ = 8 .28 x lO^T^cm 3 s" 1 (8) 

\ 9nm e J 

The force per unit volume on the neutral fluid due to ion drag is given by an expression of 
the form (Shu 1992) 

fni = ~lPPi(v - Vi), (9) 

where pi is the ion mass density, v the neutral velocity, Vi the ion velocity, and 7 is the 
so-called drag coefficient, 

= = 2 ?5 x 10l3cm 3 s -l g -l (1Q) 

rrii + p: 

Our numerical value differs somewhat from Draine et al. (1983) because we use the 
potassium value for the ion mass mj. 

Electron-neutral coupling results in a finite electrical conductivity for the disk gas 
(Krall & Trivelpiece 1973): 

2 2 

^e = = 1 r— , 11 

m e v eN m e n{av) eN 

where e > is the magnitude of the electron charge, n e is the electron number density, and 
v e N is the electron-neutral collision frequency. The associated resistivity is 



7/ 



c 



2 



4-7TOV \n, 



<c 



234(-)T 1 / 2 cmV 1 , (12) 



where c is the speed of light. The dimensionless measure of the relative importance of 
resistivity is the magnetic Reynolds number, 

Reu = (13) 
77 

where va is the magnitude of the Alfven velocity 

VA = ^ (14) 
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and B is the magnetic field vector. 

Finally, the cyclotron frequencies of the electrons, ions, and neutrals are denoted 

u) ca = (15) 

m a c 

where a is respectively e, i, and // (m^ = p). Although the neutrals are of course not 
directly affected by the magnetic field, the neutral cyclotron frequency is an important 
characteristic frequency of the analysis. 



2.2. Equations 

The fundamental dynamical equations are mass conservation, 

dp 



. )t +V-(pv) = (16) 



the dynamical equation of motion, 



dv ( B 2 \ ( B \ 

p- + (pv-V)v = -V (P + — - P V$ + --V B (17) 



and the induction equation, 

dB 



Of 



Vx (v e xB -rjVxB) 



Here, v is the velocity of the neutrals, and v e is the electron velocity, and $ is the central 
gravitational potential. 

In equation (|T?p, which is dominated by the neutral component, we have made implicit 
use of the low inertia limit for the ions. Ion momentum balance is set by balancing the 
ion-neutral drag force with the Lorentz force, 

-J x B = f ni , (19) 

c 

since all other terms in the ion momentum equation are small. It is actually f n i that 
appears on the right hand side of equation (|17); the low ion inertia approximation consists 
of instantaneously balancing this force against the Lorentz force. 



2.3. The Induction Equation 



Note the appearance of the electron velocity v e in the induction equation. In ideal 
MHD, the distinction between the fluid velocity and the electron velocity is unimportant. 
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But in fact, in the absence of resistance, it is the electron fluid in which the magnetic field 
is "frozen," since the electrons are the mobile charge carriers. The difference between ion 
and neutral velocities is responsible for ambipolar diffusion, and the difference between ion 
and electron velocities gives rise to the Hall effect (e.g., Wardle & Konigl 1993). 

In general, 

v e = v + (v e - Vi) + (vi - v) = v 1 — , (20) 

n e e 'jppiC 

We are assuming that any ions present are singly ionized (which is appropriate for a very 
low ionization gas), so that the electron and ion number densities are equal. In the ideal 
MHD limit, both v e — and — v are assumed to be negligibly small. Retaining them in 
the full induction equation leads to 

dB _ ( _ Anr,J JxB (J X B) X B\ 

— = VX i)Xfl + - '- , 21 

at \ c n e e c^pip ) 

where we have used the standard MHD approximation 

J = -^V X B. (22) 

47T 

The four terms on the right side of the induction equation can be arranged by order 
of magnitude. Reading from left to right, we denote them / (Inductive), O (Ohmic), H 
(Hall), and A (Ambipolar). Assume that typical fluid velocities are of order va, and that 
typical gradients are of order an inverse scale-height 1/h. Then, the relative sizes of these 
terms are 

1 H Lu ce A u ci 

1 Re M O Uen H 7p 

The first ratio just reflects the familiar fact that when the magnetic Reynolds number 
approaches unity, ohmic losses are important. The results of §2.1 give 

1 , o „ ,n-14 ( n \ f°S 



I Re M \n e J \v A 




4.2X10- 4 - ^ ^ ^ , (24) 



and 



H ujce^ / 8x 10 17 \ 1/2 (va\ A uj^^ / 9x 10 12 \ 1/2 / T_\ l ' 2 /v_a\ 



(25) 



Under interstellar conditions 1/ ' Rcm is tiny, even though n/n e may be large. (Note that in 
equation [13], h would be a characteristic interstellar length, not a disk thickness.) In this 
case, neither the O nor the H term is important. But the ambipolar diffusion A term may 
be non- negligible. This leads to a regime of ideal MHD in the ions, which then couple the 
magnetic field to the neutrals via ion-neutral collisions. 

In protostellar disks, our interest here, the far smaller ionization fraction can reduce 
ReM to the point where resistivity can affect the dynamics. This Reynolds number need not 
be of order unity in an MHD-turbulent fluid before it is significant; some field configurations 
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are affected when Re M is larger than 10 4 (Fleming, Stone, & Hawley 2000). The key point 
here, however, is not the well-known result that resistivity is important in protostellar disks, 
it is the less appreciated fact that HEMFs will generally be important in a magnetized disk 
whenever resistivity is. The two should be treated on the same footing (W99). 

The importance of the Hall term relative to Ohmic losses is noteworthy because it 
is surprisingly general. As formulated above, for a given ratio of magnetic to thermal 
energy, the relative importance of the Hall term is independent of the ionization level, and 
depends only upon the number density of the dominant neutrals. Furthermore, the relative 
importance of the Hall term is not a feature unique to gases of low ionization fraction. 
Using a nominal resistivity value of 5 x 10 12 T~ 3 / 2 cm 2 s _1 for a fully ionized plasma (Spitzer 
1962), we find 



o-^-^m) lirj ■ (26) 

where v e i is the electron-ion collision rate. We see that if Ohmic dissipation is of interest 
(as it would be for magnetic reconnection), the temperature and density regimes of ionized 
accretion disks imply that the Hall effect cannot be ignored. 

At typical disk protostellar disk densities, ambipolar diffusion can be neglected, though 
it may become important in the outermost regions of the disk. Since 

(J x B) x B = (J ■ B)B - B 2 J, (27) 

if the currents and fields are perpendicular, ambipolar diffusion acts like a field-dependent 
resistivity. This is the case for the special vertical field geometry treated by W99 and in 
§4.2 below: all of ambipolar diffusion may be incorporated into an effective resistivity. But 
for the more general field geometries we consider, no such simplification is possible. We 
shall accordingly assume a neutral number density much in excess of 10 13 cm -3 , and ignore 
ambipolar diffusion. The form of the induction equation we shall use henceforth is 

dB „ / „ c(V X B) x B\ 

- = Vx Ub-^xB-I^ . (28) 



3. Uniform Rotation 

Hall currents introduce novel elements into the MHD of astrophysical disks, and to 
develop an intuitive understanding, we consider a very simple problem: Hall-modified 
Alfven waves in a uniformly rotating disk threaded by a vertical magnetic field, B = Be z . 
Here we shall ignore finite resistivity and vertical structure, work in the Boussinesq limit, 
and use standard cylindrical coordinates (R, <f>, z) with the origin at the disk center. Finally, 
we consider plane wave perturbations that depend only upon z. Linearized quantities 
(indicated by 5 notation) are proportional to exp(iut — ikz), where cu is the angular 
frequency and k is the vertical wave number. 

Under these circumstances, pressure, density, vertical velocity and vertical magnetic 
field perturbations all vanish. We solve for 5v R , 5v^, 5B R and SB^. The linearized radial 
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and azimuthal equations of motion are 



ikB 

- iu5v R - 2n8v.fi - — 5B R = 0, (29) 

ikB 

- iu5vs + 2Q5v R - - — 5Bs = 0. (30) 

47rp 

The same components of the linearized induction equation are 

- iu8B R + 5Bs - ikB 5v R = 0, (31) 

Anen e 

k 2 Be 

- iu5B^ SB R - ikB 5vs = 0, (32) 

4iren e 

Note the symmetry of the problem, in particular the introduction of a "Coriolis term" into 
the magnetic field equations. These, of course, are the electromotive forces of the Hall 
effect. 

The combination k 2 Bc/47ren e will be recognized by students of plasma physics as 
the electron whistler frequency. At characteristic wavenumbers satisfying kvA ~ ^, this 
frequency is ~ Q(Q/uj c ^)(n/n e ). Hence, the ionization fraction below which the whistler 
branch couples with the disk dynamics is (n e /n) ~ (Q/uj cfJ ). 

The "whistler drift" of the field lines with respect to the bulk of the fluid is caused by 
the motion of the current-bearing electron charge carriers. Perturbing the magnetic field 
induces local currents, and it is the electron motion in the currents relative to the ions 
that causes the field lines to move relative to the ion-neutral fluid. (Recall that ambipolar 
diffusion is ignored, so the ions are locked with the neutrals.) Note the sense of the field line 
drift: if f2 and B are both upward (say), the induced whistler circular motion of the field 
lines is right-handed, following the electrons, opposite to motion induced by the dynamical 
Coriolis force. 

How does all this affect the wave response of the fluid? Consider first the case of 
aligned B and S7. The sense of the field line motion is to slow the dynamical epicycles. 
The easiest way to see this is to note that in the absence of rotation, equations ( p9]) and 
(p0|) show that the perturbed velocity vector is proportional to the perturbed magnetic field 
vector. The whistler waves would impart a circularly polarized component to the velocity 
response. This response is counter to the direction that the dynamical epicyclic motion is 
trying to impart, and hence it effectively lowers the Coriolis force. As for the magnetic 
tension, the additional Hall response is 90° out-of-phase with respect to this force, and we 
therefore expect these two to add in quadrature. If, in addition, the dynamical epicycle is 
slowed by Hall currents, the return magnetic tension force should effectively increase. If the 
field and angular velocity are counter aligned, the signs of these effects should switch. 

These expectations are born out in the dispersion relation, which is most transparently 
written 

lu 2 ± 2ficj(l - k 2 v 2 H /4tt 2 ) - k 2 {v\ + v 2 H ) = 0, (33) 
where we have introduced the Hall velocity defined by 
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Note that although this quantity has dimensions of (velocity) 2 , it may be positive or 
negative. It is positive (negative) if f2 and B are oriented in the same (opposite) sense. We 
address more general orientations in subsequent sections. 

In the dispersion relation (]3"3"|), the upper + sign corresponds to left-handed polarization 
relative to the magnetic field, the lower — sign to right-handed polarization. Whistler waves 
exist only for the right-handed branch, and correspond to the limit Q,n e — > 0. There are 
no instabilities in a uniformly rotating disk, since there is no free energy source. 

We postpone until later sections (which include differential rotation) a full discussion 
of the solutions to the dispersion relation, but this is a convenient spot to make contact 
with waves in a magnetized plasma. In the limit Q — > 0, the dispersion relation becomes 

uj 2 ± uj^-^- + k 2 v\ = 0, (35) 
4vm e e A v ' 

This is precisely the dispersion relation for a uniformly magnetized plasma with the 
displacement current and electron inertia ignored (e.g., Krall & Trivelpiece 1973). (The 
Alfven speed in the above is dominated by the neutrals, whereas in standard fully ionized 
plasma treatments it is determined by the ions, of course.) Without loss of generality, we 
may take uo > 0, allowing k to determine the propagation direction. The positive frequency 
solutions are 

k 2 Bc ( k*B 2 c 2 2 \ 1/2 
U = T ^e + \u^ + kVA ) ■ (36) 

For small wavenumbers (low frequencies) both of the above solutions reduce to Alfven 
waves. At large wavenumbers, right-handed waves (+ sign) go over to the high frequency 
whistler wave branch, whereas large k left-handed waves (— sign) are cut off at a frequency 
given by 

a; (cut off) = — (^)=M-X (37) 
fic n n 

analogous to the ion-cyclotron cut-off in a fully ionized plasma. The characteristic ionization 
fraction at which the cut-off frequency drops below the disk rotation frequency is once again 
(n e /n) ~ ft/uj cfl . 



4. Differential Rotation 



4.1. Preliminaries. 

We consider next the local stability of a differentially rotating disk threaded by a weak 
vertical field. We assume that finite resistivity and Hall currents are both present; we 
neglect ambipolar diffusion. As in the previous section, we restrict ourselves to plane wave 
disturbances of the form exp(ikz — iut). In the Boussinesq limit, this corresponds to fluid 
displacements in the plane of the disk, so vertical structure is unimportant. 

This problem has been considered by W99, but it bears re-examination for several 
reasons. It is an important baseline for understanding the effects of Hall currents in 
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protostellar disks, and is worth understanding from more than one perspective. Our 
approach emphasizes dynamical couplings rather than conductivity tensor formalism, and 
is physically quite distinct. It is also somewhat simpler. Finally, we shall examine a 
qualitatively new feature of the Hall effect that has not been discussed before: HEMFs 
can destabilize outwardly increasing differential rotation. In the process we show that the 
maximum growth rate of the instability is always the local Oort A value of the rotation (in 
magnitude), whether the angular velocity is increasing or decreasing outward. 



4.2. Axial Fields and Wave numbers 



4-2.1. Stability 

The key issue of the analysis is the local stability criterion of a disk with an axial field. 
Whereas the dispersion relation is somewhat complex, the limit u — > is simple, and this is 
relevant for understanding local stability. (In principle, local overstability may occur, but 
for self-consistency we wish to focus on nonpropagating, evanescent, local modes. For these, 
the transition between stability and instability must proceed through the point u> = 0.) 

Including the effects of differential rotation and finite resistivity along with the Hall 
terms, we find that the linearized, uo = 0, radial and azimuthal equations of motion are now 

ihB 

-2tt5vt-—5B R = 0, (38) 

*^-™SB, = 0, (39) 
where k is the epicyclic frequency, defined by 

«* = 4n> + ^. (40) 
The same components of the linearized induction equation are 

k 2 rj5B R + dB^ - ihB 5v R = 0, (41) 

47ren e 

Notice that in equation fl42|) , the whistler frequency Hall term introduces a coupling that 
is equivalent to changing the shear rate dVL/dR. It arises becauses the <\> component of 
the perturbed electron velocity (5v e ) differs from the dominant neutral velocity by a term 
involving 5J^. This produces a radial field component, 5B R . The effect is present whether 
or not there is differential rotation, but it is most striking when such motion is present, 
because it couples exactly like dQ/dR. In effect, it is the magnetic "epicyclic" term, 
combining a Coriolis-like coupling with a background angular velocity gradient. 
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If we ignore resistivity for the moment, our problem decouples into two very simple 
subproblems, one with 8v$ coupled only to SB R , and one with Svr coupled only to SB^. 
The torque must vanish, so equations (j39|) and ((4l|) combine to give 



2 2 K 2 k 2 V 2 H 

kVA+ AQ 2 



SB, = 0, 



(43) 



or the tidal and excess centrifugal forces must vanish, in which case equations 
yield 

/ dtt 2 x 



\dlnR 

If we are not to have a trivial solution, then 



+ k 2 (v 2 A + v 2 H ) SB 



R 



0. 



and (ID 

(44) 



2„,2 



Va + 



K V 



H 



An 2 



' dQ 2 
dlnR 



+ k 2 (v A + v 2 H ) =0. 



(45) 



In the standard MRI, the torque is always purely Alfvenic, and therefore it always has the 
same sign: a restoring negative torque for a positive angular displacement. This means 
that instability depends only upon the direction of the excess centrifugal force. The Hall 
effect renders the problem more interesting, allowing an interplay between the torque and 
centrifugal force. An inwardly directed centrifugal force can destabilize if accompanied by a 
positively directed torque. 

The inclusion of resistivity couples SBr and SB^, but gives only a slight modification 
to the above result: 



2„,2 



dtf 



AVL 2 ) \d\nR 



H 



+ k\v 2 A + V 2 H ))+K 2 V 2 k 2 



0. 



(46) 



This is the desired critical stability condition. To determine the sign for instability, one 
may take the limit vfj,r] —* 0, which returns us to the simple MRI. It then easily follows 
that the left hand side should be negative for instability, that is 



k 2 v\ 



;i + x)(i + 



2 

AD 2 ' 



+ 



K 2 r) 2 



< 



i + 



2 

K X 

AV? 



dQ 2 
dlnR 



(47) 



where 



x 



JL 



(4* 



is the dimensionless Hall parameter. The terms in square brackets incorporate the additional 
physics of HEMFs and resistivity; they become unity for the standard x — * MRI (Balbus 
& Hawley 1991). The dimensionless resistivity parameter may be expressed in terms of the 
magnetic Reynolds number: 

I - ^O*""' (49) 
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The behavior of the polynomial 

2 2 2 

D(x) = (l + x)(l + — ) + (50) 
determines the nature of the instability. We may write the instability criterion 

sgn{D)k 2 v 2 A < -± '—— j '- . (51) 

It is not difficult to show that if 

, 1 fQ 2 \ fv A \ rflnfi 2 

there will be a finite range of x for which D(x) < 0. In this window, all wavenumbers, no 
matter how large, are unstable, even in the presence of resistivity: the right hand side of 
the inequality (|47|) is positive, while the left hand side is negative. Clearly, the zeroes of D 
are critical for understanding the behavior of disks in the Hall regime. 

For a Keplerian disk, the window occurs when x falls approximately between —1 and 
—4 for large Reu- Roughly speaking, there is a match between the response frequency of 
the magnetic tension and the "epicyclic" frequency of the field line drift through the fluid 
caused by HEMFs. These tend to cancel each other's dynamical response. Qualitatively, it 
is as though the restoring tension were not present as a radial force, leaving the destabilizing 
dynamical tide to do its work. Within the window, going to larger wavenumbers does not 
change the sign of the radial forces on a fluid element. The azimuthal torque, however, is 
still dominated by the Alfven term, and this allows for angular momentum transfer. 

In figure (P), we show a stability plot for an r\ = Keplerian disk in the x, (kv A /fl) 2 
plane. The shaded regions correspond to instability. The instability window (i.e., the region 
in which all wavenumbers are unstable) for —4 < x < — 1 is evident, together with an 
abrupt transition for x < —4. For > x > — 1 the Hall currents destabilize by allowing for a 
larger range of destabilizing wavenumbers than would be present in the magnetorotational 
instability. For x > 0, Hall currents are a stabilizing influence. 

Figure (@) illustrates the modified instability range for large Rcm- The window is 
narrowed as shown. 

Figure (|3|) illustrates the effect of increased resistivity on stability. The window has 
now disappeared, and as Reu diminishes, an ever decreasing portion of wavenumber space 
is unstable. This region corresponds to very large wavelengths, which may exceed the global 
disk scales. This corresponds to resistive stabilization. In general, when resistivity becomes 
dominant, the instability criterion (51) leads to kv A < 0{Re~ h }). But the the right hand 



side of the inequality is a function oT~r, and at its maximum we find kv A < 0(ReM — 1/2). 
This is considerably easier to satisfy, though it corresponds to tuning the Hall parameter. 
This effect was also noted in W99. 

Figure (f|) is presented to illustrate a point of principle. It is the stability diagram 
of an 7] = 0, k 2 = 5Q 2 (i.e., dQ 2 /dR > 0) disk. There is an instability window, analogous 
to its Keplerian counterpart, for — AQ 2 /k 2 > x > — 1, abrupt stability for x > —4Q 2 /k 2 , 
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and instability for an ever diminishing wavenumber domain as x falls below —1. In the 
presence of sufficiently large Hall currents, with counteraligned axial angular velocity and 
magnetic field vectors, the instability criterion is dQ 2 /dr ^ 0. Any differential rotation law 
is potentially unstable. 



4-2.2. Dispersion Relation 



We next obtain the full dispersion relation. Assume that the time dependence of the 
perturbations takes the form exp(crt). (This form keeps the coefficients of the dispersion 
relation real.) The linearized radial and azimuthal equations of motion are now 



a5vR — 2£l8vri 



ikB 

.. . ikB 
a5v, + -8v R - — 



SB f 



0. 



SBj, = 0. 



The same components of the linearized induction equation are 



{a + k 2 n)5B R + 



k 2 Bc 



[a + k 2 n)5B 4 



47ren L . 
k 2 Bc dVL 



SB^ — ikB 8v R = 0, 

5B R — ikB Svs = 0, 



\4:iien e dlnR 
The resulting dispersion relation is 

= a 4 + 2r]k 2 a 3 + C 2 a 2 + 2r]k 2 (K 2 + k 2 v 2 A )a + C t 
where the constants C 2 and C are given by 



k + 2k 2 v\ + r] 2 k A + 



k 2 v 2 H , dn 2 



Ail 2 K d\nR 



+ k 2 v 2 H ) 



and 



Co 



2 , k2v H 

Va + An 2 



dQ 2 
dlnR 



+ k 2 (v 2 A + v 2 H ) 



+ n 2 n 2 k A 



(53) 
(54) 

(55) 
(56) 

(57) 

(58) 
(59) 



Balbus & Hawley (1992a) conjectured that the maximum growth rate of any 
instability feeding off the differential rotation in a disk is given by the local Oort A value, 
o A = (l/2)\dQ/dlnR\. All magnetic field configurations of the magnetorotational instability 
have this value for their most rapid growth rate. However, the A value conjecture went 
beyond Alfven tension as the destabilizing agent, it suggested that whatever the proximate 
cause, linear perturbations can grow no faster than a a- The reasons are rooted in the 
dynamics of the differential rotation process itself, not in magnetism per se. Now, it is far 
from obvious that in the absence of resistivity, the rather unwieldy dispersion relation ( |57D 
has precisely this value as its maximum growth rate. But it does, as we now demonstrate. 
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To begin, write equation 
We define 

s = a/fl, k = k/Q, X = {kvA/Q) 2 

Note that 



in dimensionless form, with all rates normalized to Q. 

Y = {kv H /Q) 2 . (60) 



K 



4 + 



dlntt 2 



dlnR 



The f] = dispersion relation is 



s 4 + 



~ 2 nx, Y fdlnn 2 ^ 
k 2 + 2X + - — — + Y 
4 \ d\nR 



s 2 + [X + 



Yk 2 \ (d\ntt 2 



d\nR 



+ X + Y 



(61) 



(62) 



At the maximum growth rate s = s m , partial differentiation of the above with respect to X 
and Y gives the two equations 



X + \(k 2 + A) 



j(K 2 + 4) + j(n 2 + s 2 m 



2 ldlnQ 2 
' Sm ~2~d~^R 



(63) 



lrflnfi 2 
4 d\nR 1 



~2 , ,2 
K, -+- S m 



(64) 



Using the identity (^) and eliminating Y between equations (|63|) and (|64D , leads after 
simplification to a remarkable result, 



(X + S 2 m + K< 



°m 



1 /dlnfl' 
16 V rflni? 



2\ 2 



(65) 



Since the first factor must be positive definite, we conclude that 



1 

16 



d\nR 



(66) 



This is equivalent to 



i.e., the Oort A value! 



dlnQ 



d\nR 



(67) 



We must now verify that this value of s 2 ^ is a solution of the dispersion relation (|62|) 
for X and Y satisfying the system of equations (|63D and (|64"D. It is easily shown that with 
the Oort A value chosen for s m , the determinant of this system vanishes. We thus use only 
one of these equations, which we take to be (|63D. First we substitute for s 2 ^, and solve for 
X in terms of Y. This leads to 



X = — 



Idlntt 2 Y 



- —\K 



2 dlnR 8 



4) 



1 /dlnQ 2 
16 V d\nR 



(68) 



Note that in the Y = limit of this equation, X > is possible only for decreasing outward 
angular velocity profiles. The magnetorotational instability requires dfl 2 /dR < 0. The 
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addition of the Hall effect allows for the possibility that Y < 0, and thus for a perfectly 
well-defined positive X parameter, even if dfl 2 /dR > 0. That both senses of angular 
velocity gradient may be treated on the same footing, accords with the well-known result 
that the state of minimum energy for a disk of fixed angular momentum is solid body 
rotation (Lynden-Bell & Pringle 1974). In principle any angular velocity gradient could be 
unstable, since lower energy states consistent with angular momentum conservation exist. 
A dynamical path to these low energy states is required however, and that is what HEMFs 
can provide. 

Replacing X by the above expression wherever it occurs in equation (|62|) yields, after 
some simplification, 

s 4 + c 2 s 2 + c = 0, (69) 

with 



1 



C-2 



(Y - 4) s 



1 fd\nQ 2 \ 2 



8 V dlnR 



and 



c 



1 fdlntt 2 ^ 
64 V d\nR 



1 fdlntt 2 * 
4 \ d\nR 



(Y-Af 



The standard solution is 



'd\nSF 
d\nR 



{Y - 4) 2 ± (Y - Af 



(70) 



(71) 



(72) 



Choosing the + sign returns the Oort A solution, exactly as desired and consistent with 
the Balbus & Hawley (1992a) conjecture. (W99 obtained the Keplerian upper limit by 
numerical solution of the dispersion relation, and noted that it was the same as the standard 
magnetorotational instability.) 



Graphical solutions in the XY plane are shown in figures ([51) and 



4.3. General Axisymmetric Disturbances 

We now consider the axisymmetric behavior of the instability with more general field 
geometries and wavenumbers. We shall ignore buoyancy, so our analysis holds either for 
a strictly polytropic disk, or locally at the midplane for any constituative relation. The 
perturbation wavevector is 

k = k R e R + k z e z , (73) 

and disturbances have space-time dependence exp(?fc ■ r + at). The linearized equations 
are now that of mass conservation (Boussinesq limit): 

k R 5v R + k z 5v z = 0, (74) 

and the equations of motion, 

a5v R - 2Q6V4 - 6B R + % ^-{B^ 5B+ + B Z 5B Z ) + ik R — = 0, (75) 
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*6v* + ^5v R - '-^p- 6B+ = 0, (76) 

ik R B R ik z m x „ , n \ , **z*-P n / 77 n 

otfvs j + - — [Bfj, oB^ + B R oB R ) H = 0. (77) 

47rp 47rp p 

If one considers a polytropic disk, then 5P/ p may be replaced directly by a pure enthalpy 
perturbation, SQ, say. Otherwise it is to be taken at face value, a pressure perturbation 
divided by the density p. Ultimately, either choice leads to the same result, because the 
term serves only as a place holder. 

The linearized induction equations are 

(a + k 2 V )5B R + < k - B "> kz 5B _ . B)8v R = 0, (78) 
Airen e 

(a + k 2 r])5B z - ^ ' B ^ R 5Bs - i(k • B)8v z = 0. (80) 
4:7ien e 

The dispersion relation that emerges after a straightforward but lengthy effort is 

a 4 + 2 V k 2 a 3 + C 2 a 2 + 2 V k 2 (^k 2 + (k • v A f J a + C = 0, (81) 



with 



rb O/luiO/ig \ (X 111 A Zi/lo/tg z/ 



and 



(83) 

A sufficient criterion for instability is clearly 

C < 0. (84) 

This would ensure that the sign of the quartic on the left side of equation (|8lf) changes as 
a passes from small to large positive values. Hence there would have to be a positive root 
somewhere. The necessity of (|84"D as an instability criterion is more difficult to prove, but is 
very likely to be true. It is straightforwardly proven in the limit of vanishing resistivity that 
Cf > leads to stability (since the analysis reduces to the roots of a quadratic), and the 
presence or resistivity will almost certainly cause further stabilization. Numerical analyses 
are all consistent with (|8~4" ) as a necessary and sufficient criterion for instability. 

The most important difference between strictly axial geometry and the present 
calculation lies with the form of the coupling of the Hall term and wavenumber. Whereas 
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the axial result makes it appear that the Hall term stabilizes or destabilizes according 
to whether f2 • B is positive or negative, equations (|81|)-(B3|) show that the coupling 
is actually (k ■ f2)(fc • B). The significance of this form ofthe coupling is that under 
conditions of marginal stability, if there is any finite radial field component, there will 
always be wavenumbers which make this term negative (destabilizing). This is of particular 
importance for the nonlinear development of the instability, where power is cascaded 
throughout the wavenumber spectrum. 

Consider the case of astrophysical interest, dQ 2 /dR < 0. If the final factor of equation 
(jS3|) (the torque term) is negative, the preceding factor (the centrifugal term) will also be 
negative, and obviously no instability will be possible. Let us assume, therefore, that the 
torque is positive. Then, a necessary condition for instability is 

which differs from the standard MRI only by an additive Hall term on the left hand side. 
(This is a necessary condition because it omits resistivity.) The effect of this is shown 
schematically in figure (7). We have introduced the dimensionless Hall parameter, 

Ha = cMfc ' f g) , (86) 

the characteristic whistler frequency normalized to Q. 

Let us finally note that HEMFs will tend to stabilize geometries with vanishing axial 
field. In the absence of the Hall term, the standard Bz = MRI, has its maximum growth 
rate for disturbances with kz — > oo. But with Hall EMFS included, arbitrarily large axial 
wavenumbers will not allow Cq < 0. Instead, these modes become stable waves. 



4.4. Nonaxisymmetric Disturbances 



4-4-1- Shearing Coordinates 

One of the most important reasons for obtaining this general axisymmetric formula 
(jSll) is that it provides a basis for understanding the behavior of local nonaxisymmetric 
disturbances, to which we now turn. 

We shall work in the local shearing coordinate system introduced by Goldreich & 
Lynden-Bell (1965) in their pioneering study of local gravitational disk instabilities, and 
used in a magnetic context by Balbus & Hawley (1992b) and Terquem & Papaloizou (1996). 
The shearing coordinates are defined by 

R' = R, <j)' = <j) - Q(R)t, z' = z, (87) 

which leads to 

d d dVt d d d d d , , 

-t- 



dR dR' dRdch' dd> ddf dZ dZ> 
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The partial time derivative must be taken holding shear coordinates constant, which is just 
the Lagrangian time derivative for the undisturbed disk flow: 

Local disturbances depend on these sheared coordinates as exp[i(k' R R' + mcf)' + k z z')], where 
the wavevector k' is a constant. This means that we may effect the transformation from 
Eulerian to (unperturbed) Lagrangian coordinates simply by replacing the static Eulerian 
wavevector k R with 

k R (t) = k' R -mt^, (90) 

and calculating the time evolution with the Lagrangian derivative. Throughout this section, 
k shall denote the shearing wavevector (kn(t), m/R, kz)- 

In the presence of shear, the toroidal magnetic field grows linearly with time: 

B^t) = B^(0)+tB R ^. (91) 
Note, however, that k • B is time-independent: 

k • B = k' R B R + + k z B z . (92) 

This has important consequences for nonaxisymmetric behavior. 



4-4 -2. Linear Equations 
The linearized dynamical equations are 

Tfl 

knSv R + —Sv^ + k z Sv z = 0, (93) 
K 

d$VR ., (8P B-8B\ .fk-B\. n 

d5v z ., f5P BSB\ .{k>B\ rn 

alt 2Q K R \ p Att P J \ 4vrp J K ' 

For the induction equations, it is helpful to define an auxiliary vector A by 

A = -r]k 2 SB + C(fc ' B \ kx SB). (97) 
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The magnetic field equations are then 

d5B R 
dt 

d5B z 
dt 



i(k> B)5v R + A R , (98) 



i(k- B)5v z + A z , (99) 



<U)D " ^{k^B)5v <t> + -^-5B R + A <t> . (100) 



dt v ' v dlnR 
The above equations imply that 

= -,* 2 (fc • S), (101) 

so that if chosen to vanish initially, k • SB will vanish throughout the evolution in a 
numerically stable manner. If k • SB vanishes, then k • A — 0. 

Because of the time dependence of the wavevector, the nonaxisymmetric problem is 
more complex than its axisymmetric counterpart. The approach we shall use is to reduce 
the system to two coupled second order differential equations for 5B R and SB^, which are 
then solved numerically (Balbus & Hawley 1992b). This is an algebraically tedious process, 
and we present only the final result: 



d_ (d5B L _ A ' 
It \ dt r , 



+ 



2mk r dVL ( d5B r 



k 2 dr 



dt 



2rQ k 2 z 
m k 2 

-aA - 



2mQk z fd5B 



rk 2 \ dt 



A.) 



'd5B z 
dt 



-A, 



(k • v A f SB r , 



(102) 



d (d5B z 



dt \ dt 



-A, 



m k 2 



(d5B z \ 
k z \ dt A j+k r 



^2mk z d(rQ) ( dSB r ^ 



rk 2 dr 



dt 



(d5B r " 



(k-v A ) 2 8B z . 



(103) 



These equations are formally the same as equations (2.19) and (2.20) of Balbus & Hawley 
(1992b), but with the first-order time derivatives of 5B transformed to 



d5B d5B 



dt 



dt 



A. 
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4-4-3. Behavior of Perturbations 



Although equations ( |102| ) and ( |103| ) appear rather opaque, it is possible to understand 



the onset of instability by focusing on the comparatively simple axisymmetric inequality 
flSSp. In essence, we find that the nonaxisymmetric evolution unfolds as a series of 
axisymmetric problems. In going from axisymmetry to nonaxisymmetry, the only term in 
the Co stability criterion that acquires a time-dependence is k. Consider the evolution of a 
strongly leading disturbance, /cr large and negative. As kn initially moves towards zero, the 
factor (kz/k) 2 goes from much less than, to nearly unity (assuming m 2 /(kzR) 2 is small). 
The right hand side of ([85]) starts out very small, and grows to its maximum value of 
\dVt 2 / d\nR\. At this stage, the formal instability criterion must be satisfied, or it will never 
be satisfied over the entire course of the evolution. Then, as kn grows from zero to large 
and positive, the right hand side of ( p5] ) once again diminishes with time, and eventually 
the flow will stabilize. 

The situation is summarized in figure (|7|), which shows the stability and instability 
regions of a Keplerian disk in the (k/kz), (k • va/Q) 2 plane. We have adopted the 
nominal values Ha = —1.6 and Ha = 0.4, and for comparison have also displayed the 
magnetorotational instability results (Ha = 0) of Balbus & Hawley (1992b). Note that 
despite the complexity of the full system of equations, in the end the HEMF simply slides 
the hyperbolic stability curve up and down along the ordinate. Since k • va is constant, 
the evolving path of a perturbation is along a horizontal line in this plane, since only 
the abscissa coordinate changes with time. The nature of the unstable response is then 
determined entirely by how long this horizontal trajectory remains inside the instability 
region. 

In figure (||]), both the Hall effect and resistivity are included. From equations (|83| ) 
and (|84]), we see that the marginal stability curves in the (k/kz), (k • «a/^) 2 plane depend 
on the resistivity only through the parameter r]k 2 z /VL. We have adopted r]k 2 z /Vt = 1. To 
orient oneself, this would correspond to a disk with aspect ratio h/R = 0.1, k z R = 100, 
and (cs/va)Rgm = 100. We see clearly in figure (|8|) that the resistivity stabilizes small 
radial wavelengths. In figure (^), to allow for direct comparison, we have plotted both the 
resistive and nonresistive cases. 



Figu res (J7|)~(|9|) have been obtained strictly by numerically solving equations ( |102j ) 
and ( |103| ). Nevertheless, each of these curves is almost indistinguishable from those 
obtained simply by setting C = 0, where C is defined by (1531). Therefore, although the 



criterion (|84|) has been derived for axisymmetric perturbations, it is very accurate for 



nonaxisymmetric disturbances also, with k/kz varying with time. Note that the curves 
displayed in figures (0)-© depend on m only implicitly, through the ratio k/kz- The role 
of m is limited to determining the evolution of k/kz, not the properties of the stability 
curves. 

In figure (flOl) , we show the evolution of the perturbed radial field 5Br (the curves 
for 5Bz and SB^ are similar) for the same Ha values as above, and for both r] = and 
rjk'z/Q = 1. The value of (k • va) 2 /Q 2 = 1-5 is chosen because it lies in the unstable regime 
for an axisymmetric 77 = disk, for all the Hall parameters. Initially k/kz is very nearly 
unity, and the perturbation grows until the wavevector moves out of the relevant unstable 
region in the (k/kz), (k • va/Q) 2 plane. When 77 = 0, explosive growth is achieved for 
the three values of Ha considered. As expected, the cases Ha = —1.6 and Ha = 0.4 are 
respectively more and less unstable than the case Ha = 0. The effect of the resistivity is to 
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suppress the growth of the perturbation at large wavenumbers, corresponding to very early 
and very late times in the evolution. The dramatic reduction of peak amplitudes is evident. 
The Ha < case, however, remains quite unstable even with r]k z /Vl = 1. 

To conclude, we have found that if m/kzR is small, the nonaxisymmetric behavior of 
linear perturbations may be understood on the basis of axisymmetric behavior. (Large m 
disturbances are stable.) Whether growth occurs or not depends upon the instantaneous 
location of the wavenumber in the (k/kz), (k • Va/Q) 2 plane. Shear causes the wavenumber 
to retrace a horizontal line in this plane, first from right to left (for an initially leading 
distrubance), then the reverse. If this path takes the wavenumber into an axisymmetrically 
unstable zone, the disturbance grows. The effect of resistivity, not surprisingly, is to dampen 
growth, with the smallest wavelengths being the most affected. 



5. Summary 

In a gas which is at least partially ionized, the magnetic field is, but for resistivity, 
frozen into the electron fluid. The difference in the mean electron velocity and the 
center-of-mass fluid velocity gives rise to HEMFs in the gas. Somewhat surprisingly, in 
astrophysical environments these can often be more important than Ohmic dissipation, 
even in a fully ionized plasma. This point is often not appreciated. In protostellar disks, the 
Hall effect is small compared with resistivity only if the neutral density is much in excess of 
10 18 cm 3 , or if va "C c§. 

The dynamical consequence of HEMFs is the appearance of whistler waves as a mode 
of the gas response. These waves are carried nominally by the electrons, but inductive 
electron-ion coupling and ion-neutral collisional coupling together assure that the bulk of 
the fluid is involved. It is the interplay between these right-handed circularly polarized 
whistlers and the Coriolis driven epicycles that affects the stability of magnetized disks. 
Whistler waves become dynamically important when the ionization fraction drops below 

The Hall effect allows disks with either decreasing outward or increasing outward 
angular velocity profiles to become unstable. By way of contrast, the standard MRI affects 
only those disks with a decreasing outward profile, even though any angular velocity 
gradient increases the rotation energy for fixed angular momentum. The maximum growth 
rate is always found to be the Oort A value, (1 /2)\dQ/ din R\. 

Whether HEMFs stabilize or destabilize a disk depends not upon fl • B, but upon 
(k • Q)(k • B). The latter should be negative for destabilization. If B has a radial 
component, then it is always possible to find a wavenumber to make this parameter 
whatever value is desired. Thus, Hall physics will always tend to make some wavenumbers 
more unstable. Whether or not the disk actually become unstable will depend on the value 
of the resistivity. 

The nonlinear consequences of HEMFs are likely to be particularly important for our 
understanding of transitions between active and quiescent states in both dwarf novae and 
FU Orionis outbursts. Neither of these systems is presently understood at a fundamental 
level. Numerical simulations which include only Ohmic losses have yet to demonstrate that 
MHD turbulence can be sustained at the ionization fractions that may be present in these 
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cool disks (Fleming et al. 2000, Menou 2000). Numerical simulations including both the 
Hall effect and Ohmic dissipation should prove most informative. 
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Fig. 1. — Stable and unstable regions in the (Icva/^) 2 , x = (v h /va) 2 parameter space for the 
case of zero resistivity. The black diamond corresponds to the location of the critical axial 
wavenumber for the MRI. Positive values of v 2 H always stabilize because the right-handed 
polarization of whistler waves decreases the effective magnitude of the destabilizing shear. 
Negative values at first destabilize by effectively increasing the shear rate, but ultimately 
stabilize in the form of whistler waves at large magnitudes of the Hall parameter. Notice 
that for —4 < x < — 1, all wavenumbers are unstable; magnetic tension is eliminated as a 
stabilizing agent. 
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Fig. 2. — As in fig. (1), but with finite resistivity, K 2 rf/v\ = 0.35. The window of instability 
at all wave numbers is still present even with resistance, but it is more narrow. 
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Fig. 3. — As in fig. (2), but with larger resistivity, K 2 rf jv\ = 1. The window of instability at 
all wave numbers has vanished, and only the largest wavelengths remain formally unstable. 
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Fig. 4. — Instability diagram of an rj = 0, k 2 = 5f2 2 disk profile, corresponding to an 
outwardly increasing rotation profile. No instability is present in such a disk without the 
Hall EMF. With it, the disk is unstable, and shares what seems to be the universal maximum 
growth rate, 0.5d£l/dlnR. 
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Fig. 5. — Growth rate in the X, Y plane for an 77 = Keplerian disk. (See text for definitions.) 
Only regions of instability are shown, with the height proportional to the growth rate. The 
vertical axis has been suppressed, but the maximum growth rate of the "ridge" is 0.75 fl. 
Y < corresponds to counter- aligned S7 and B. 
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Fig. 6. 
profile. 



. — As in fig. (5), but with k 2 = 5Q 2 , a disk with an outwardly increasing velocity 
Only Y < regions can be unstable. 
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Fig. 7. — Stability and instability regions in the (k • va) 2 /VL 2 — (k/kz) plane for 
nonaxisymmetric disturbances, which correspond as well to axisymmetry. The Hall 
parameters Ha are: Ha = (solid line), Ha = —1.6 (dotted line) and Ha = 0.4 (dashed 
line). For Ha = and Ha = 0.4, the unstable region is under the curve. For Ha = —1.6, it 
is between the two curves, k • va is constant for a shearing wavevector. For strongly leading 
disturbances, k/kz is initially large, and the point corresponding to the wavevector moves 
to the left in the plane on a constant k • Va line. The minimum value of k/kz is unity (to 
order m 2 /(k z R 2 )). After attaining its minimum, the wavector point retraces its path to the 
right. For values of k • va smaller than some critical value which depends on Ha, a finite 
portion of time is spent in the unstable region, and substantial growth may occur. 
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Fig. 8. — Same as in Fig. |7| but with a non zero resistivity r]k 2 z /Vt = 1. Here, for each 
value of Ha, the unstable region is between the curve and the vertical axis. The resistivity 
stabilizes the large k,R disturbances. 
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Fig. 9. — Same as in Fig. [7] and |8] but here, to allow for comparison, we have plotted both the 
cases rj = (solid line) and r]k 2 z /Vt = 1 (dotted line) on the same panel. The Hall parameters 
are as shown. 
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Fig. 10. — Evolution of the perturbed radial field SB R for a Keplerian disk with 
(k ■ va) 2 /^ 2 = 1.5. The initial amplitudes are 5B R = 0.3 and 5B Z = SB^ = 0, and the 
wavenumbers are k z R = 100, m — 1. The initial value of k R R is determined by V • SB = 0. 
The left panels have zero resitivity, the right panels have r\k 2 /Q = 1. Hall parameters are as 
shown. 



